An Extension of the Complete Model-Based 
Decomposition of Polarimetric SAR data 


Abstract—The model-based decomposition originated from 
Freeman-Durden three-component decomposition (FDD) has 
been widely applied in polarimetric synthetic aperture radar 
(POLSAR) data processing for its clear physical interpretation 
and easy implementation. Numerous improvements have been 
proposed to settle the two main drawbacks of FDD, i.e. incomplete 
utilization of the polarimetric information in the coherency 
matrix and the negative scattering power problem. Recently, Cui 
et al. proposed a complete model-based three-component 
decomposition of POLSAR coherency matrix data which has 
successfully settled the two aforementioned drawbacks. However, 
the three scattering components’ powers are not totally derived 
using scattering models, and the remained coherency matrix 
(RCM) obtained by subtracting the volume scattering component 
from the coherency matrix is not consistent with the models of 
surface scattering and double-bounce scattering components. As 
an extension of Cui’s method, this letter is dedicated to develop a 
novel method to discriminate the surface scattering and 
double-scattering components from the RCM both using 
scattering models. With the orientation angle (OA) variation and 
helix angle variation compensated for the RCM, the RCM is 
automatically consistent with the models of surface scattering and 
double-scattering components. The OA variation and helix angle 
variation compensation for the RCM is done by unitary 
transformations of the eigenvectors of the RCM. Then the 
compensated RCM is used to derive the powers of surface 
scattering and double-scattering components using scattering 
models. The powers of surface scattering and double-bounce 
scattering components are positive. The effectiveness of the 
proposed method is demonstrated by processing POLSAR data. 


Index Terms—¥reeman-Durden three-component 
decomposition, helix angle variation, nonnegative eigenvalue 
decomposition, orientation angle variation, polarimetric synthetic 
aperture radar. 


I. INTRODUCTION 


Polarimetric target decomposition plays an important role in 
target identification, target classification, and geophysical 
parameters retrieval. Huynen pioneered the work of target 
decomposition. He decomposed the distributed target into a 
single target and the target “noise” [1]. The instability, 
non-uniqueness, as well as the preference for symmetry and 
regularity, made that the Huynen decomposition has not been 
widely used [2]. Nowadays the hotspots of target 
decomposition research mainly concentrate on two categories, 
i.e. the Cloude-Pottier eigenvalue/eigenvector decomposition 
and the model-based decomposition originated from 
Freeman-Durden three-component decomposition (FDD). For 
the Cloude-Pottier eigenvalue/eigenvector decomposition, the 
polarimetric coherency matrix is decomposed into three 
orthogonal eigenvectors which represent three orthogonal 
single targets, respectively [3], while for the FDD, the 


covariance matrix is decomposed into three scattering 
components, i.e. surface scattering component, double-bounce 
scattering component, and volume scattering component [4]. 
We only focus on model-based decompositions in this letter. 
Because reflection symmetry (RS) is assumed, FDD only 
accounts for five of nine independent parameters of the 
coherency matrix, i.e. the T(1,3) and T(2,3) of the coherency 


matrix T are not used, and thus polarimetric information loss 
is unavoidable. Due to blindly assigning all the 
cross-polarization to the volume scattering component, the 
negative power of surface scattering or double-bounce 
scattering component may occur in FDD. In fact, numerous 
improvements [5-16] have been achieved to solve these two 
problems. Yamaguchi et al. [5] added a helix component to 
account for the off-diagonal elements of the coherency matrix, 
and six independent parameters are accounted for. The negative 
power problem is also alleviated in some degree as the helix 
component also accounts for some cross-polarization power. 
An et al. [6] and Lee et al. [7] introduced deorientation 
operation to the coherency matrix to mitigate the negative 
power problem, and reduced the independent parameters from 
nine to eight. The Yamaguchi’s four-component decomposition 
with deorientation operation can account for six out of eight 
parameters [8]. Based on RS assumption, nonnegative 
eigenvalue decomposition (NNED) was proposed by van Zyl et 
al. [10, 11] to avoid the negative power problem, but the 
polarimetric information is still not fully utilized. 

Recently, Cui et al. [12] proposed a complete model-based 
three-component decomposition. The volume scattering 
component is determined by the NNED without RS assumption, 
and thus full utilization of polarimetric information is 
guaranteed. In [12], two algorithms are introduced for 
discriminating the second and third components, i.e. surface 
scattering and double-bounce scattering components. In the 
first algorithm, the powers of these two components are derived 
as the eigenvalues of the remained coherency matrix (RCM) 
after interpreting the scattering mechanisms of the de-oriented 
eigenvectors of the RCM. However, the powers are not derived 
using scattering models as mentioned by Lee et al. [14]. In the 
second algorithm, the RCM is first de-oriented, and from which 
a maximum single scatterer, surface scatterer or double-bounce 
scatterer, is derived using scattering model as the second 
scattering component. But, the third component is not derived 
using scattering model, and the derivation procedure is 
complex and ambiguous in physical meaning as mentioned by 
An et al. [15]. An et al. [15, 16] considered that the (3, 3) 
element of the coherency matrix of the third scattering 
component should be zero, because for a pure surface scatterer 
or double-bounce scatterer the (3, 3) element is zero which can 
be verified by the models of surface scattering and 


double-bounce scattering components: 
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Thus, An et al. [15] improved Cui’s method by implementing 
deorientation to minimize the (3, 3) element of the coherency 
matrix of the third scattering component, and the (3, 3) element 
can be much closer to zero. However, the (3, 3) element is still 
not zero which is still not consistent with the models of (1) and 
(2). 

As an extension of Cui’s method, this letter is intended to 
develop a novel method to discriminate surface scattering and 
double-bounce scattering components from the RCM, and the 
proposed method manifests two main advantages: 1) all the 
three components are derived all using scattering models; 2) 
After compensating the orientation angle (OA) variation (two 
OAs) and helix angle variation (two helix angles) for the RCM, 
the RCM is consistent with the models of (1) and (2), which 
leads to that all the powers are concentrated on co-polarization 
elements. Only in this circumstance, the compensated RCM is 
suitable to be used to derive the powers of the surface scattering 
and the double-bounce scattering components using the 
scattering models of (1) and (2). The powers are all positive as 
demonstrated by the results obtained by the proposed method 
using POLSAR data. 

We briefly review the Cui’s method in Section II, introduce 
the proposed method in Section II, present experimental 
results and some discussions in Section IV and Section V, and 
finally conclude this letter in Section VI. 


II. BRIEF REVIEW OF CUI’ S METHOD 


The complete model-based three-component decomposition 
proposed by Cui et al. can be represented as [12] 


which can be also rewritten as 
T —- yTy =R; + Polp. (4) 


All the matrices in (4) should be physically meaningful, i.e. 
the matrices should be positive semi-definite and do not have 
negative eigenvalues. The matrix T—P,Ty should also be 
positive semi-definite. Both Cui et al. and Wang et al. showed 
that the maximum possible volume scattering power is the 
minimum root of the cubic equation [12, 13] 
where det(-) denotes matrix determinant. The method of 
volume scattering power derivation is based on NNED without 
RS assumption. The derived maximum possible volume 
scattering power has also been suggested as the reasonable 


value to be used [10, 11]. 


The RCM, 7’ =T —TyP, , contains two single scatterers, i.e. 
the surface scatterer and the double-bounce scatterer. Cui et al. 
proposed two methods to discriminate these two single 
scatterers, which will be briefly reviewed in the following. 

A. The Algorithm 1 

The first algorithm is based on eigenvalue/eigenvector 

decomposition of the RCM. The RCM is decomposed as 

T' = Ak! + Akk}! (6) 
where the superscript H denotes the complex conjugate 
transpose. The eigenvectors kų and k are then de-oriented, 


and are finally discriminated as surface scattering mechanism 
or double-bounce scattering mechanism based on the 
co-polarized phase. The corresponding eigenvalues 7, and A, 


will be assigned as the powers of the corresponding scattering 
mechanisms. 
B. The Algorithm 2 


The second algorithm tries to fit a maximum power single 
surface (or double-bounce) scatterer with azimuth slope 
modulation to the RCM T’. The RCM can be represented as 


T' = RR” (O)TsR(Q)+T" (7) 
where R? (OTsR(0) is the azimuth slope modulated single 


scatterer, 0 is the OA, and the OA rotation matrix is 
1 0 0 
R(@)=|0 cos20 sin28 |. (8) 
0 -sin2@ cos20 


The value of 0 can be derived by maximizing Ry . The RCM is 
OA rotated as 

R(O)T'R" (0) = PTs + R(O)T'R" (0) (9) 
The scattering mechanism of the power maximized single 
scatterer is determined by the value of the parameter @ of Ty. 


If | B | <1, a power maximized rotated surface scatterer is fitted 


to the RCM, and a power maximized rotated double-bounce 
scatterer is fitted when | B | <1. The last scattering component is 


determined using the coherency matrix R(O)T'R” (A) , whose 


corresponding scattering matrix is calculated, then de-oriented, 
and at last the scattering mechanism is determined based on 
co-polarized phase. 

It can be seen that the powers of surface and double-bounce 
scatterers are derived from eigenvalues and eigenvectors for 
Algorithm 1, not from scattering models. For Algorithm 2, 
although the second scattering component is derived using 
scattering model, it is not the case for the third scattering 
component derivation, and the derivation procedure is 
complicated and ambiguous in physical meaning [15]. 


IHI. THE PROPOSED METHOD 


As an extension of Cui’s method, the proposed method is 
devoted to discriminating the surface scattering and 


double-bounce scattering components from the RCM. 

First, the volume scattering power is derived using Cui’s 
method, and then the RCM is decomposed as shown in (6). 

Then, the OA variation and helix angle variation are 
compensated for the RCM, i.e. the two eigenvectors are 
individually compensated for its OA and helix angle. The 
eigenvectors are modeled by the Kennaugh-Huynen’s 
co-diagonalization of the scattering matrix [17, 18]. S; is 
assumed to be the corresponding scattering matrix of 
eigenvector k;, i=1,2. 

Si = Rs(G)Rs (Tj) Sai Rs (Ti Rs (-G), i = 1,2 (10) 

where Sy; is the diagonal matrix whose diagonal elements are 
complex co-eigenvalues x; and x; , respectively. 0, and T; 


are the OA and helix angle, respectively. The corresponding 
—sin 6; 
cos 6; 


l , respectively. Then we 


cos 6; 
unitary transformation matrices are Rg (6) =| . 
sin 6; 
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expand and rewrite (10) in Pauli basis, and consequently, the 
corresponding eigenvector k; can be represented as 
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Hence, the OA 0; can be easily derived from (11) as 
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The OA compensation for k; is carried out by 
where R(6;) is the OA rotation matrix of (8), and k; is OA 
compensated eigenvector as 
T 
ee a Pes A nay SE as 
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The helix angle 7; is easily derived from (14) as 


ik; (3 
Ti = lian J il ) : 
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The helix angle compensation for eigenvector is carried out as 
kj =U _k; (16) 


where k;' is the OA and helix angle compensated eigenvector. 


(12) 


(15) 


The helix angle unitary transformation matrix is 
cos2r 0 jsin2r 

U, = 0 1 0 The OA and helix angle 
jsin2t 0  cos2r 


compensated eigenvector has the form 
T 
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We can see that the cross-polarization element of the 
eigenvector is zero now after OA and helix angle 
compensations. 

The compensated RCM T; can be synthesized by combining 


(17) 


the two compensated eigenvectors together as 


2 
Ti =X Akik}. (18) 
i=l 


It is interesting to see that not only T¿(3,3) but 
also Ty (1,3), T: (3,1), T2 (2,3), T: (3,2) of the compensated RCM 
are all zero now. It is noteworthy that T% is consistent with (1) 
and (2), and T% can be used to derive the surface and 


double-bounce scattering components’ powers using scattering 
models of (1) and (2) 
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are thus determined. The proposed extension of Cui’s method is 
presented as Algorithm 1 below. 
Algorithm 1: An Extension of the Cui’s Method 

: Input T and Ty 

: Solve cubic equation: det(T-PyTy)=0 

: Select the minimum root to be the value of Py 

: Obtain the RCM: T' =T - yTy 


: fori=1to2do 
: Obtain OA 6; with (12) 
Implement OA compensation for k; with (13) 
Obtain helix angle 7; with (15) 
10: Implement helix angle compensation for k; with (16) 


1 
2 
3 
4 
5: Implement eigendecomposition: T' = kk" + Akak} 
6 
7 
8 
9 


11: Obtain OA and helix angle compensated k;' 
12: end for 


13: Obtain compensated RCM: Ty = Akik” + Apk5k4” 
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16: end if 
17: Output: Ps, Pp, Py 


IV. EXPERIMENTAL RESULTS 


In this section, the decomposition results of Cui’s Algorithm 
1 and the proposed method on RADARSAT-2 dataset of San 
Francisco area are presented and compared. As the 
performances of Cui’s Algorithms 1 and 2 are similar [12], the 
decomposition results of Cui’s Algorithm 2 are omitted here for 
saving space. 

The used RADARSAT-2 C-band data of San Francisco area 
were acquired on April 9, 2008 in a fine-beam and 
quad-polarization mode. The original data are in single-look 
complex scattering matrix format, and the azimuth and ground 
range resolutions are about 4.82 m and 4.73 m, respectively. 
The data were then 6-look processed both in azimuth and 
ground range directions to get the coherency matrix data, after 
which they have a dimension of 2402 X 470 pixels. There are 
ocean, park, and built-up areas in the scene. Color-coded 
decomposition results of Cui’s Algorithm 1 and the proposed 
decomposition are shown in Fig. 1. We only present a portion 
of the whole scene. The surface scattering, double-bounce 
scattering, and volume scattering contributions are color coded 
in blue, red, and green, respectively. It can be seen from Fig. 1 
that, the color of ocean areas, which are surface scattering 
dominant areas, is blue as expected. The color of park area, 
which is volume scattering dominant area, is green. The color 
of built-up areas, which are double-bounce scattering dominant 
areas, is red in general. But in some largely oriented built-up 
areas, the color is yellow or even green, which means that the 
volume scattering contribution is still over-estimated even with 
the proposed method and Cui’s method implemented. The 
average span normalized 7'(3,3) and average span 


normalized 7: (3,3) by the proposed method for the whole 


image are computed, and are 6.59% and zero on the whole 
image, respectively. As pointed out in [15, 16], the T'(3,3) of 
the RCM, i.e. the cross-polarization power of surface scattering 
and double-bounce scattering components represented by (1) 
and (2), respectively, should be zero, which has been achieved 
by the proposed method. 

To specifically demonstrate the different characteristics 
between the proposed method and Cui’s Algorithm 1 further, 
three typical patches are selected for detailed processing. The 
selected three patches are ocean area, built-up area, and central 
park area as shown in Fig. 1(b). The surface scattering, 
double-bounce scattering, and volume scattering power 
proportions are shown in Table I. As can be seen from Table I, 
the power proportions by the two methods are similar generally, 
because of which that the visual effects of the decomposition 
results by the two methods of Fig. | are similar generally. As 
can be seen from Table I, the surface scattering, double-bounce 
scattering, and volume scattering contributions are dominant in 
ocean area for patch 1, in built-up area for patch 2, and in park 
area for patch3, respectively. One difference between the 
proposed method and Cui’s Algorithm 1 is that the proposed 
method gets a little higher surface scattering power, which is 
because that further helix angle compensation is adopted by the 
proposed method. Comparing (14) and (17), we can find that 
the cross-polarization element becomes zero after further helix 


angle compensation implemented, and the corresponding 
power of the cross-polarization will be transferred from 
T'(3,3) to T'(1,1) of the coherency matrix. The T'(,D 
represents the surface scattering contribution in a sense. Thus, 
the surface scattering power obtained from the proposed 
method is a little higher. 

In reality, it is still a problem in verifying that the result of 
which method is consistent with the ground truth. Hence, we 
just compare the performances of the proposed method and 
Cui’s Algorithm 1 here to demonstrate the characteristics of the 
two methods and the effectiveness of the proposed method. 


V. DISCUSSION 


A. The Necessity of Helix Angle Variation Compensation 


In [1], Huynen used one parameter F to represent the helicity 
of a distributed target, and the helicity is kind of global shape 
twist of the target. In [5], Yamaguchi et al. introduced helix 
scattering component as the fourth component of the scattering 
contribution from a target. In [17, 18], Huynen and Touzi both 
used helix angle to characterize the asymmetry of the target. 
Thus, the helicity is the inherent characteristic of the target. 

When using scattering models to discriminate the surface 
scattering and double-bounce scattering components from the 
RCM, the models of (1) and (2) are not consistent with the 
RCM [15, 16]. However, if we use the proposed method with 
OA variation compensation and further helix angle variation 
compensation for the RCM, then the compensated RCM is just 
consistent with the models (1) and (2). Actually, except for 
OAs of the surface scatterer and the double-bounce scatterer, 
the helix angles of the surface scatterer and the double-bounce 
scatterer also contribute to the cross-polarization, which can be 
demonstrated by and the third element of (14) and the helix 
scattering component of the Yamaguchi four-component 
decomposition [5] 


0 0 0 
Thelix = 0 1 +1 (20) 
0 F1 1 


When (1) and (2) were being modeled in [4], any OA and helix 
angle were not considered. Thus, OA variation and helix angle 
variation of the RCM must be compensated for before it is used 
to discriminate the surface scattering and double-bounce 
scattering components. 

In addition, the helix angle 7} and 7 of the whole image is 


generally small. Thus, the effect on the decomposition results 
of helix angle compensation of (16) is not so obvious, but the 
proposed method has necessarily accommodated the RCM to 
the scattering models. 


B. Interpretation of the Two Eigenvectors 


Questions may be raised that how to interpret the two 
eigenvectors and which one is of surface or double-bounce 
scattering mechanism. It is noteworthy that we do not interpret 
the scattering mechanisms of the two eigenvectors. The two 
orthogonal eigenvector can be considered as the projections of 
all the surface scatterers and double-bounce scatterers in one 


resolution cell to the orthogonal axes expanded by the two 
orthogonal eigenvectors. We just compensate for the OA 
variation and helix angle variation via Cloude’s 
eigenvalue/eigenvector decomposition of the RCM, and 
combine the two compensated eigenvectors together as (18) to 
get the compensated RCM which is used to derive the powers 
of surface and double-bounce scattering components. The 
compensated RCM is also a positive semi-definite Hermitian 
matrix and physical realizable. 


C. Future Work 


As shown in the results, although with the proposed method 
and Cui’s method implemented, the scattering mechanism is 
still misinterpreted in some largely oriented built-up areas. The 
reason may be that the used volume scattering models do not fit 
all the scattering cases. The adaptive volume scattering model 
fitting all the scattering situations needs further study. 


VI. CONCLUSION 


An extension of the Cui’s complete model-based 
decomposition is proposed in this letter, which provides a new 
method to discriminate the surface and double-bounce 
scattering components via compensating the OA and helix 
angle variation for the RCM using scattering models. As an 
alternate to Cui’s method, it accommodates the RCM to the 
models of surface scattering and double-bounce scattering 
components. Its effectiveness is verified on the real POLSAR 
data. 
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TABLE I 
AVERAGE SCATTERING PROPORTIONS OF SURFACE SCATTERING, 
DOUBLE-BOUNCE SCATTERING, AND VOLUME SCATTERING 
COMPONENTS IN SELECTED THREE PATCHES 


Algorithm Ps/SPAN Pp/SPAN _ Py/SPAN 
Patch 1 Proposed method 94.38% 2.52% 3.1% 
Cui’s Algorithm 1 94.35% 2.55% 3.1% 
Patch 2 Proposed method 31.49% 57.49% 11.02% 
Cui’s Algorithm 1 29.19% 59.79% 11.02% 
Patch 3 Proposed method 36.66% 15.23% 48.11% 
Cui’s Algorithm 1 35.57% 16.32% 48.11% 


Fig. 1. Decomposition results on RADARSAT-2 data of San Francisco by (a) 
Cui’s Algorithm 1 and (b) the proposed method. 


